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Abstract. We assess transport properties of heavy quarks in the Quark-Gluon Plasma (QGP) using static 
heavy-quark (HQ) potentials from lattice-QCD calculations in a Brueckner many-body T-matrix approach 
to evaluate elastic heavy-quark-light-quark scattering amplitudes. In the attractive meson and diquark 
channels resonance states are formed for temperatures up to ~1.5T C , increasing pertinent drag and diffu- 
sion coefficients for heavy-quark rescattering in the QGP beyond the expectations from perturbative-QCD 
calculations. We use these transport coefficients, complemented with perturbative elastic HQ gluon scat- 
tering, in a relativistic Langevin simulation to obtain HQ pt distributions and elliptic flow (V2) under 
conditions relevant for the hot and dense medium created in ultrarelativistic heavy-ion collisions. The 
heavy quarks are hadronized to open-charm and -bottom mesons within a combined quark-coalescence 
fragmentation scheme. The resulting single-electron spectra from their semileptonic decays are confronted 
with recent data on "non-photonic electrons" in 200 AGeV Au-Au collisions at the Relativistic Heavy-Ion 
Collider (RHIC). 



1 Introduction 

One of the most interesting questions in high-energy nu- 
clear physics is that about the properties of the hot and 
dense medium created in ultra-relativistic heavy-ion col- 
lisions. Finite-temperature lattice-QCD (1QCD) calcula- 
tions of strongly-interacting matter predict a phase transi- 
tion from hadronic matter to a quark-gluon plasma (QGP) 
at a critical temperature, T c ~ 180 MeV [T]. In the re- 
cent years the experimental program at the Relativistic 
Heavy-Ion collider has resulted in convincing evidence for 
the formation of such a hot and dense partonic state [HO 

HP]. 

The heavy charm and bottom quarks are particularly 
valuable probes for the properties of this medium since 
they are created in the primordial hard collisions of the 
nuclcons within the colliding nuclei. Thus, they form a 
rather well defined initial state and interact with the hot 
and dense fireball during its entire evolution. Recently, 
measurements of the transverse-momentum distributions 
of "non-photonic single electrons" (e^), which originate 
mainly from the semi-leptonic decays of open-charm and 
-bottom mesons, in 200 AGcV Au-Au collisions at RHIC 
have found a surprisingly large suppression at high trans- 
verse momenta (pt) (i.e., a small nuclear modification fac- 
tor, Raa) and a large elliptic-flow parameter, v%. Both 
findings indicate that during the lifetime of the hot and 
dense fireball heavy quarks come close to thermal equilib- 
rium with the medium [SHU]. 



The theoretical challenge is to understand the corre- 
sponding thcrmalization times of heavy quarks from the 
underlying microscopic scattering processes with the con- 
stituents of the QGP, in particular how the heavy quarks, 
despite their large masses, vtlq ^> T c , become part of 
the collective flow of the fireball. In calculations of the 
pertinent transport coefficients from perturbative QCD 
(pQCD), based on gluon-brcmsstrahlung energy loss, in- 
cluding elastic HQ scattering, one has to artificially tune 
the coupling strength beyond the applicability range of 
perturbation theory [91110]. It has also been shown that 
the convergence of the perturbative series for the HQ dif- 
fusion coefficient is quite poor [11] . Thus, non-perturba- 
tive approaches have to be used to explain the strong HQ 
couplings necessary. One suggested mechanism is the for- 
mation of D- and i3-meson resonance excitations in the 
deconfined phase of QCD matter [T2lll3] which has lead 
to a quite satisfactory description of the data at RHIC. 

This paper is organized as follows: In Sec. [2] we use 
HQ static potentials from lattice-QCD calculations at fi- 
nite temperature in a many-body Brueckner T-matrix ap- 
proach to calculate elastic HQ light-quark scattering-ma- 
trix elements in the medium [T411T5] . We show that after 
inclusion of a complete set of color channels, taking into 
account 1 — 0,1 states in the partial-wave expansion of 
the T-matrix, the resonance states, conjectured in the ear- 
lier approaches, are confirmed by these interactions, which 
are in principle free of tunable parameters. The resulting 
clastic-scattering amplitudes are used in Scc.[3]to calculate 
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drag and diffusion coefficients for a Fokkcr-Planck equa- 
tion [T6 l [T2 l [T7] . describing the rescattering of the heavy 
quarks within the hot and dense sQGP fireball. In the next 
step we employ a relativistic Langevin simulation to find 
the corresponding HQ p t distributions, using a thermal- 
fireball parameterization, including elliptic flow for non- 
central heavy-ion collisions. To confront these spectra with 
the e ± data from the PHENIX and STAR collaborations 
at RHIC, in Sec. [4] we use a combined quark-coalescence 
and fragmentation model to hadronize the heavy quarks to 
D and B mesons which then are decayed semi-lcptonically 
leading to the final spectra which can be directly con- 
fronted with recent data on nonphotonic single electrons 
in 200 AGcV Au-Au collisions at RHIC. The paper closes 
with brief conclusions and an outlook (Sec. [5J. 

2 HQ scattering in the QGP 

In this Sec. we calculate in-medium matrix elements for 
elastic scattering of heavy quarks (Q = c, b) with light 
quarks q = u, d, s in a Brueckner-like many-body ap- 
proach, assuming that a static heavy-quark light-quark 
potential, V{r), can be employed as the interaction kernel. 
Such a model has been used in the vacuum to successfully 
describe D- meson spectra and decays [I8jll9j. Further, we 
assume that the effective in-medium potential can be ex- 
tracted from finite-temperature 1QCD calculations of the 
color-singlet free energy F\(r,T) [SOpT] for a static QQ 
pair as the internal potential energy by the usual thermo- 
dynamic relation [TH221I23U24] . 



U 1 (r,T) = F 1 (r,T)-T 



diq(r, T) 
df ' 



(1) 



For application as a scattering kernel in a T-matrix equa- 
tion, the potential as to vanish for r — > oo. Thus we choose 
the accordingly subtracted internal potential energy, 



Vi(r,T) = C/i(r, T) — Ui(t ► oo,T). 



(2) 



In 1QCD simulations one finds that U\(r — ► oo,T) is a 
decreasing function with temperature which could be as- 
sociated as a contribution to the in-medium HQ mass, 
m Q{T) = tuq + U\(r — > oo,T)/2 where mo denotes the 
bare mass. However, this leads to problems since close to 
T c the asymptotic value, U±(r — > oo,T), develops a pro- 
nounced peak structure. Thus, in this calculation, we as- 
sume constant effective HQ masses, m c = 1.5 GeV and 
m b = 4.5 GeV. 

We also consider the complete set of color channels for 
the Qq (singlet and octet) and Qq (anti-triplet and sextet) 
systems, using Casimir scaling as in leading-order pQCD, 



i. 



1, 



(3) 



which is also justified by recent 1QCD calculations of the 
finite-T HQ free energy [551125] . 

This approach is in principle parameter free in the 
choice of the interactions, since their strength is taken 
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Fig. 1. (Color online) Diagrammatic representation of the 
Brueckner many-body scheme for the coupled system of the 
T-matrix based on the 1QCD static internal potential energy 
as the interaction kernel and the HQ self-energy. 



from first-principle 1QCD simulations. However, there arc 
considerable uncertainties in the potentials (a) between 
different lattice calculations and (b) in the extraction and 
parameterization of the corresponding free energies, par- 
ticularly their temperature dependence needed to subtract 
the entropy term in Eq. (jTJ) . In addition, the very notion 
of an "in- medium potential" is not a unique concept [27], 
and its identification with the internal potential energy 
may be seen as an upper limit in interaction strength. We 
use three different parameterizations of Fx 2-! 2211X71 : 

[Wo] of quenched 1QCD [20], 
[SZ] of two-flavor 1QCD [28], and 
[MR] of three-flavor 1QCD [59]. 

The resulting potentials from [Wo] and [SZ] are compara- 
ble to a numerical extraction from three-flavor 1QCD [53] , 
while that from [MR] is deeper than the other two for 
T < 1.6T C , but falls off faster at higher temperatures. 
The resulting uncertainty in the transport coefficients (see 
Sec. [3]) amounts to up to 40%. To define the Brueckner- 
type many-body scheme the four-dimensional (4D) Bethe- 
Salpetcr (BS) ladder approximation, symbolized in dia- 
grammatical form by the upper panel of Fig. [JJ has to 
be reduced to a 3D Lippmann-Schwinger (LS) equation, 
neglecting antiparticle components in the quark propaga- 
tors, in order to implement the static potential from 1QCD 
via Eqs. ([JJ3]). After this reduction the LS equation in the 
color channel, a £ {1, 3, 6, 8} reads [i"4pl 

T a (E;q',q) =V a (q',q)-J A±_V q (q' ,k)G qQ (E;k) 



(2tt) 3 

T a {E-k,q)[l-f F {Lol)-f F {^)] 



(4) 



with the Fourier-transformed potentials, 

V a (q\ q)= I d 3 rK(r) exp[i( Q - q')r]. (5) 

Further, E, q and q 1 denote the energy and incoming and 
outgoing momenta in the center-of-mass (CM) frame, re- 

1 Here and in the following all vertex and Green's func- 
tions are understood as the retarded real-time quantities which 
can be derived as analytic continuations of the correspond- 
ing imaginary-time (Matsubara) quantities of thermal quan- 
tum field theory. 
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spcctivcly. 



Jf = 



1 



(6) 



exp(uj/T) + 1 

is the Fermi-Dirac distribution. The quark-dispersion re- 
lations are determined in quasi-particlc approximation by 



k 2 



(7) 



where for simplification we do not solve the fully self- 
consistent scheme in Fig. Q] but use a fixed mass of m q = 
0.25 GeV, m c = 1.5 GeV, and m b = 4.5 GeV for the 
light, charm, and bottom quarks, respectively. Finally, the 
two-particle-gQ propagator in ([¥]) is given in terms of the 
Thompson-reduction scheme |30j 



G q Q(E; k) 



1 



1 



(8) 



with a quasi-particle width for both light and heavy quarks 
of -2Sf Q = 0.2 GeV. 

The solution of the LS equation (01 is simplified by 
using a partial-wave expansion of the potential and T- 
matrix, 

V a {q', q) = 4n ^(2Z + l)V a ,itf , q)Pi[cos Z(q, q% 
I 

T a (E; q', q) = 4tt^(2Z + l)T a>l (E; q' , g)P z [cos Z(q, q% 



(9) 



which leads to the ID LS equations, 



T atl (E;q',q) = V a ,i(q',q) + - / dkk 2 V a ,i(q' , k)G Qq (E; k) 



xT ajl {E-k,q)[l-f F {u2)-f F {ul% 

(10) 

for the partial- wave components, T a j, of the T-matrix, 
which are solved numerically with the matrix-inversion al- 
gorithm of Haftel and Tabakin [31] . We restrict ourselves 
to S (I = 0) and P (I = 1) waves. As can be seen from 
Fig. [21 in the dominating attractive color-singlet Qq and 
color-antitriplct Qq channels, close to the critical tempera- 
ture, T c ~ 180 MeV, resonance states above the threshold, 
Ethr = "iQ + m q are formed, similar as conjectured in [121 
1 1 3| - However, in this full in- medium scheme the resonances 
melt at higher temperatures T > 1.7T C and T > 1.4T C , re- 
spectively. As we shall sec in the next section, contrary to 
the expectation from perturbative calculations, using the 
[Wo] parameterization of the potential, this leads even to 
decreasing transport coefficients with increasing tempera- 
ture, i.e., the decreasing interaction strength of the poten- 
tial overcompensates the higher density of the medium. 

The HQ self-energy, diagrammatically represented by 
the lower panel of Fig. [I] is then given by 



ds\d a f k 2 dkdx 



6 J 4tt 2 
x T a (E;p,k), 



[/-F(wfc) + /s(w + a> fe )] 



(11) 



sing, T=1.2 T c 
sing, T=1.5 T c 
trip, T=1.2T C 
trip, T=1.5T C 




Fig. 2. (Color online) Imaginary part of the S- wave in-medium 
T matrix for cq and cq scattering in the color-singlet and - 
antitriplet channels based on the parameterization of the 1QCD 
potential energy by [Wo]. 
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Fig. 3. (Color online) Real and imaginary parts of the c-quark 
self-energy as a function of three-momentum at different tem- 
peratures. 

where dsi = 4(2/ + 1)^V/ denotes the spin-isospin de- 
generacy of the light quarks (which is already implicitly 
assumed in the light-heavy quark interaction of our ap- 
proach, which is in line with the free D-meson spectrum 
32] ) in the Z th partial wave and d a the color degeneracy 
in the corresponding channel. Here we assume an effective 
number of light-quark flavors, Nf = 2.5, to account for 
the smaller strange-quark density. For the charm quarks, 
cf. Fig. [3] the calculation leads to an in-medium width of 
r c = — 2ImZ' c ~ 200 MeV, which justifies our simplifying 
assumption in the T-matrix calculation above. 



3 HQ transport in the QGP 

To evaluate the motion of the heavy quarks in the hot 
and dense fireball, consisting of light quarks and gluons, 
we employ a Langevin simulation of the Fokker-Planck 
equation, 



dfq 
dt 



_d_ 

dp, 



(Pilfo) 



d 2 
PiPj 



{Bijfo)- 



(12) 
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The drag or friction coefficient, 7, and diffusion coeffi- 
cients, 



By -B - p2 



P 



(13) 



are calculated from the invariant scattering-matrix ele- 
ments [16]. Taking into account elastic scattering of the 
heavy quark with a light quark or antiquark the latter 
given in terms of the above calculated T-matrix by 



J2 \ M \ 2 = ^r( s - m l + m g) 2 - m Q + m \) 

xN f J2 d a (\T a ,i=o{s)\ 2 + 3|T 0i ,=i(«) cos(# cm )| 2 ) 



2 ) 2 



(14) 



We define the averaging operator 

d 3 q 



(X( P )) 



1 

2E n 



(2^f2E q 
d 3 p> 



3„i 



(15) 



{2n) 3 2E p , lc 
(2Tr)W(p + q-p' -q')f q (q)X(p>) 



for a heavy-quark observable, X, over the elastic scatter- 
ings per unit time of the heavy quark with momentum p 
with a light quark of momentum q, changing their mo- 
menta to p 1 and q' . Here, f q is the (thermal) distribution 
of the light quarks in the medium. Then we can calculate 
the transport coefficients as 



7(|P|) = (1) 
Bo(\p\) = - 



_ {P-Pl 

<p' 2 )- <(P ' P ' )2) 



p- 



Bi(\ P \) = - 



{{P-P'?) 



P 



-2(p' ■p)+p 2 (l) 



(16) 
(17) 

(18) 



However, it turns out that, in order to guarantee the pro- 
per equilibrium limit of the heavy quarks with the me- 
dium, we have to enforce Einstein's fluctuation-dissipation 
relation for the longitudinal diffusion coefficient, which is 
used here in its relativistic form [T71I13I. 



B 1 = ET-y. 



(19) 



The nonperturbative HQ light-quark scattering-matrix ele- 
ments are supplemented by the corresponding perturba- 
tive elastic HQ gluon-scattering ones [33] . The i-channcl 
singularity is regulated by a gluon-Debye screening mass 
of m g = gT with a strong coupling constant g = y/47ta s , 
using a s = 0.4. As shown in Fig. [3J close to T c the equili- 
bration times of r cq = I/7 ~ 7 fm/c for charm quarks are 
a factor of ~ 4 larger than the values from a corresponding 
pQCD calculation, reminiscent to the results based on the 
model, assuming the survival of ZTmeson like resonance 
states above T c [T2lll3j . In contrast to this and other cal- 
culations of the HQ transport coefficients, here the drag 
coefficients decrease with increasing temperature because 
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Fig. 4. (Color online) The drag coefficient, 7, as a function of 
HQ momentum, calculated via (|16[) with scattering-matrix el- 
ements from the non-perturbative T-matrix calculation (using 
the parameterization of the 1QCD internal potential energies 
by [Wo]) compared to a LO perturbative calculation based on 
matrix elements from 1331. 



of the "melting" of the dynamically generated resonances 
at increasing temperatures due to the diminishing inter- 
action strength from the 1QCD potentials. 

To solve the Fokker-Planck equation (fT2| under con- 
ditions of the sQGP medium produced in heavy-ion colli- 
sions, we use an iscntropically expanding thermal fireball 
model, assuming an ideal-gas equation of state of Nf = 2.5 
effective massless light-quark flavors and gluons. The to- 
tal entropy is fixed by particle multiplicities at chemical 
freeze-out which we assume to occur at the critical tem- 
perature, T c = 180 MeV. For semi-central collisions the 
fireball is chosen to be of elliptic-cylindrical shape with 
isobars given by confocal ellipses with a perpendicular 
radial-flow field, scaling linearly with the distance from 
the center as seen in hydrodynamic calculations |34j to 
which also the (average) radial flow velocity and clliptic- 
ity, i>2, is fixed. To compare to "minimum-bias" data on 
e ± spectra in 200 AGeV Au-Au reactions at RHIC we 
simulate collisions with an impact parameter of b = 7 fm, 
implying an initial spatial eccentricity of about 0.6. Us- 
ing a QGP-formation time of 0.33 fm/c leads to an initial 
temperature of 340 MeV. The evolution stops after the 
fireball has undergone a mixed QGP-hadronic phase after 
about 5 fm/c, at which the radial velocity of the fireball 
has reached about a radial velocity of about v±_ = 0.5c at 
the surface and an ellipticity of V2 = 5.5%. 

Given this description of the medium, (|12p is solved 
with help of an equivalent relativistic Langevin simulation 
which is defined by the stochastic equation of motion for 
a heavy quark at position, x, and momentum p: 

Sx = ^St, 6p = --f{t 7 p)p6t + 5W{t 7 p + 5p), (20) 
E 

where 5W is a stochastic force distributed normally, 



P(SW) oc exp 



(B- 1 ) jk SW j 5W k 



(21) 
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Fig. 5. (Color Online) Nuclear modification factor, Raa, for 
central (upper panel) and elliptic flow, V2, for semicentral 
(lower panel) 200 j4GeV Au-Au collisions for charm quarks 
from the Langevin simulation, using the T-matrix results 
for the transport coefficients employing the parameterizations 
[Wo] and [SZ] for the 1QCD potentials, compared to a calcu- 
lation based on pQCD and the resonance-model interactions 
in Ml 



B 1 denoting the inverse of the diffusion-coefficient ma- 
trix (UHl). Note that in and ([21]) the diffusion coeffi- 
cients are to be evaluated at the updated momenta p + 
Sp. This Hanggi-Klimontovich realization of the stochas- 
tic process together with the dissipation-fluctuation rela- 
tion (|19p ensures the correct equilibrium limit in the long- 
time regime t ^> I/7 [35] , After evaluation of the time 
step (|20|) the resulting momenta are Lorentz boosted to 
the laboratory frame. 

The initial condition for ([20]) is given by the phase- 
space distribution of the heavy quarks. The spatial dis- 
tribution is determined with a Glauber model for heavy- 
quark production. The initial pt spectrum is determined 
from data on p-p and d-Au collisions at RH1C as follows: 
The c-quark spectra are taken from a modified PYTHIA 
calculation to fit D and D* spectra in d-Au collisions [55] , 
assuming (5-function fragmentation. After decaying this 
spectrum to single they saturate corresponding data 
from p-p- and d-Au [51157] collisions up to p t — 3.5 GeV. 
The missing yield at higher p t is assumed to be filled with 
the corresponding contributions from B mesons, leading 
to a cross-section ratio of a b i/a C c — 5- 10 -3 and a crossing 
of the c- and 6-decay electron spectra at pt — 5 GeV. 



1.8 
1.6 
1.4 
1.2 

3 1 
* 0.8 

0.6 

0.4 

0.2 





T-matrix [Wo] 

T-matrix [SZ] 

Reso, r=0.4-0.75 GeV 

■-^JS^ pQCD, <x s =0.4 



Au-Au Vs=200 GeV (central) 
bottom quarks 



2 3 
Pt (GeV) 



Au-Au Vs=200 GeV (b=7 fm) T-matrix [Wo] 
bottom quarks T-matrix [SZ] 

Reso, r=0.4-0.75 GeV 
pQCD, a s =0.4 




2 3 
P T (GeV) 

Fig. 6. (Color online) The same as Fig. [5] but for bottom 
quarks 

In Figs. [5] and [5] we show the nuclear modification fac- 
tor, Raa, defined by R AA = P Q (t Rn ,p t )/ P Q (0,p t ) (where 
tfi n denotes the time at the end of the mixed QGP-hadronic 
phase) in central and the elliptic flow, 



"2 = ({pi -pI)/pI), 



(22) 



in semicentral 200 AGeV Au-Au collisions. We compare 
the results from the T-matrix model with the parameteri- 
zations by [Wo] and [SZ] of the 1QCD potentials with the 
those using pQCD or the resonance-model interactions of 
Rcf. [T3]. While for charm quarks for the [Wo] potential 
the result for R AA is comparable to the upper end of the 
uncertainty band of the resonance-model calculation, the 
V2 is slightly enhanced at low pt- The reason for this be- 
havior is the decrease of the transport coefficients with 
increasing T: While the suppression of the pt spectra at 
high p t is due to the evolution along the whole history 
of the fireball, leading to comparable effects at the end 
of the mixed phase, the anisotropic flow is mostly devel- 
oped at the later stages and thus can be transferred to the 
heavy quarks at the end of the evolution efficiently, when 
the drag coefficient become larger due to the dynamical 
formation of the resonance states. The T-matrix result 
with the somewhat less attractive [SZ] potential leads to 
the usual ordering of the coefficients (increasing with in- 
creasing temperature) and thus shows weaker effects for 
both the R AA and v-i than the result of the resonance 
model. For b quarks the T-matrix calculations yield larger 
medium modifications of the pt spectra than the resonance 
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model which is due to the mass effect, leading to stronger 
binding effects for the resonances in the T-matrix calcula- 
tion. As to be expected, the effects of pQCD-based trans- 
port coefficients on the HQ spectra for both charm and 
bottom quarks is much weaker than the non-perturbative 
ones via the resonance-scattering mechanism. 



4 Single-electron observables at RHIC 



The last step toward a comparison of the above described 
model for HQ diffusion in the QGP with the single-electron 
Pt data from RHIC is the hadronization of the HQ spec- 
tra to D- and S-mesons and their subsequent semileptonic 
decay to e^. Here we use the quark-coalescence model de- 
scribed in [3E, 39 . In the recent years, the coalescence of 
quarks in the hot and dense medium created in heavy-ion 
collisions has been shown to provide a successful hadroni- 
zation mechanism to explain phenomena such as the scal- 
ing of hadronic elliptic-flow parameters, V2 with the num- 
ber of constituent quarks, V2,h(j>t) = nhV2 : q(pt/nh), where 
rih = 2(3) for mesons (hadrons) denotes the number of 
constituent quarks contained in the hadron, h, and the 
large p/n ratio in Au-Au compared to p-p collisions [381 
40,41]. Quark coalescence is most efficient in the low-p t 
regime where most c and b quarks combine into D and B 
mesons, respectively. To conserve the total HQ number, 
we assume that the remaining heavy quarks hadronize via 
((5-function) fragmentation. 

As shown in Fig. [7] the Langevin simulation of the HQ 
diffusion based on transport coefficients from the 1QCD 
static potentials, followed by the combined quark-coale- 
scence fragmentation description of hadronization to D 
and B mesons and their subsequent semileptonic decay, 
successfully accounts simultaneously for both the nuclear 
modification factor, Raa-, and the elliptic flow, i>2, of sin- 
gle electrons in 200 AGeV Au-Au collisions 0[7] at RHIC. 
The uncertainty due to the two different parameteriza- 
tions of the potentials by [Wo] and [SZ] is not so large. 
However, the deviations from other parameterizations are 
bigger and will be demonstrated elsewhere [42]. The ef- 
fects from the "momentum kick" of the light quarks in 
quark coalescence, an enhancement of both, Raa and V2 , 
is important for the quite good agreement of both observ- 
ables with the data. As can be seen from the lower panel 
in Fig. [3 within our model the effects from the mixing 
of the £?-meson decay contribution to the spectra be- 
comes visible in the region of p t ~ 2.5-3 GeV. A closer in- 
spection of the time evolution of the pt spectra shows that 
the suppression of high-p t heavy quarks occurs mostly in 
the beginning of the time evolution, while the V2 is built 
up later at temperatures close to T c which is to be ex- 
pected since the V2 of the bulk medium is fully developed 
at later stages only. This effect is also pronounced for the 
[Wo] parameterization of the HQ potential since in this 
case due to resonance formation the transport coefficients 
become largest close to T c . 
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Fig. 7. (Color online) Upper panel: single-electron spectra 
from the T-matrix calculation of HQ diffusion in the QGP 
based on the [Wo] (solid line) and the [SZ] (dash-double-dotted 
line) parameterizations of the 1QCD static HQ potential in 
comparison to data from the PHENIX [7] and STAR [7] collab- 
orations in 200 AGeV Au-Au collisions at RHIC. The dashed 
line shows the result when only 5-function fragmentation is 
considered for hadronization. Lower panel: Raa and t>2, as in 
the upper panel for the [Wo] parameterization for the elec- 
trons from D (dashed line) and B mesons (dash-dotted line) 
separately. 
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5 Conclusions and outlook 

We have used static potentials from finite-temperature 
1QCD calculations within a Brueckner-type many-body 
calculation, complemented by pQCD HQ-gluon elastic- 
scattering matrix elements, to assess drag and diffusion 
coefficients for c and b quarks in the QGP in a princi- 
pally parameter free approach, however plagued with large 
uncertainties in the determination of the relevant poten- 
tial from lattice data. The diffusion of heavy quarks in 
the QGP is calculated with a Langevin simulation. The 
medium is parameterized as an expanding thermal fire- 
ball (including anisotropic flow for semiccntral heavy-ion 
collisions) with an equation of state of a massless gas 
of gluons and Nf = 2.5 light-quark flavors. To confront 
this model with data on non-photonic single electrons in 
200 AGeV Au-Au collisions at RHIC, we have used a com- 
bined coalescence-fragmentation model to hadronize the 
heavy quarks to D and B mesons which subsequently de- 
cay semi-leptonically. The resulting p t spectra agree with 
recent data on the nuclear modification factor, Raa and 
elliptic flow, V2 quite well. 

In a schematic estimate from the evaluated drag and 
diffusion coefficients leading to this results, based on ki- 
netic theory using either pQCD for a weakly coupled plas- 
ma [4"3"ll4"4"] or the strong-coupling limit applying AdS/CFT 
correspondence [45"ll46] , we find values for the space-diffu- 
sion coefficient 2irTD s = 4-6 and a viscosity to entropy- 
density ratio of r\j s = 2-5/ (47r) (to be compared to the 
conjectured AdS/CFT bound (ri/s) min = 1/(4tt)), indi- 
cating a strongly coupled (liquid like) quark-gluon plasma 
close to the phase transition [47] . 

In future works detailed studies of the uncertainties in 
the potential approach is necessary, in particular about 
the question whether a static-potential approach is justi- 
fied and which in-medium potential (i.e., free or internal 
potential energy or combinations thereof) should be used 
to describe the interactions of heavy quarks within the 
sQGP. First steps in this direction have been made in [27] 
in an effective-field theory approach for non-relativistic in- 
medium QCD-bound states. Also an inclusion of inelastic 
processes like gluon bremsstrahlung which should become 
effective at higher p t is mandatory for a complete picture 
of the in- medium behavior of heavy quarks [48] . 

Another step, which is quite natural, given that within 
our model the underlying microscopic effect of the strong 
coupling of the heavy quarks to the medium is the forma- 
tion of resonance states close to T Cl is to substitute the 
quark-coalescence model with a transport-model based 
resonance-recombination model for hadronization in the 
QGP [HI] which obeys the conservation laws for energy 
and momentum as well as the second law of thermodynam- 
ics. This model shares with quark-coalescence hadroniza- 
tion the phenomenologically successful feature of the scal- 
ing of V2 with the hadrons' number of constituent quarks. 
Additionally it leads to scaling with the transverse kinetic 
energy well, provided the parton distribution and flow is 
calculated with a dynamically consistent scheme as the 
Langevin simulation discussed in the present paper |50j . 
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